Nuclear enhancement and suppression of diffractive structure functions at high 

energies 
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We compute diffractive structure functions for both protons and nuclei in the framework of Color 
Glass Condensate models with impact parameter dependence. These models have previously been 
shown to provide good agreement with inclusive F2 measurements and exclusive vector meson mea- 
surements at HERA. For nuclei, they provide good (parameter free) agreement with the inclusive 
F2 data. We demonstrate good agreement of our computations with HERA measurements on in- 
clusive diffraction. We extend our analysis to nuclei and predict the pattern of enhancement and 
suppression of the diffractive structures functions that can be measured at an Electron Ion Collider. 
We discuss how the impact parameter dependence crucially affects our analysis, in particular for 
large invariant masses at fixed Q 2 . 

PACS numbers: 13.60.Hb,24.85.+p 



I. INTRODUCTION 



The discovery that about 15% of Deeply Inelastic Scat- 
tering (DIS) events at HERA are diffractive events [J, Q 
has focused attention on the nature of hard diffrac- 
tive scattering in QCD at collider energies. In partic- 
ular, it is observed that the ratio of diffractive to in- 
clusive cross-sections is nearly independent of the en- 
ergy [!, IK, and that both cross-sections display geo- 
metric scaling d, d, 0]- It was noted some time ago 
that these features could be understood H in a sim- 
ple model, the Golec-Biernat-Wusthoff (GBW) dipole 
model [9(, which incorporated the physics of QCD sat- 
uration at high energies. The saturation of parton densi- 
ties is due to non-linear multi-parton effects such as re- 
combination and screening which deplete the gluon den- 
sity at small x [l(|. These non-linear effects are large 
for modes in the hadronic wavefunctions with transverse 
momenta k± < Q Sl where Q s (x), appropriately called 
the saturation scale, is a scale generated by the multi- 
parton dynamics. Though the Golec-Biernat-Wusthoff 
model explains qualitative features of the inclusive and 
diffractive data, it fails in detailed comparisons to the 
data. This is primarily because, except for the quark 
masses, the model does not contain geometric scaling vi- 
olations, such as for instance the bremsstrahlung limit of 
perturbative QCD (pQCD) that applies to small dipoles 
of transverse size r <C l/Q s (x). 

The GBW model is significantly improved by including 
the appropriate DGLAP behavior for dipoles with small 
transverse sizes [ll|. This "DGLAP improved" dipole 
model arises naturally in the classical effective theory of 
the Color Glass Condensate (CGC) [HG3E1. For in ~ 
stance, in the CGC [HI, [n|, one obtains the well known 
expression [r7l.[l8| for the inclusive virtual photon hadron 



cross section 
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dz 
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where 
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represents the probability for a 



virtual photon to produce a quark-anti-quark pair of size 

r = \yt\ and (ry, x, by) denotes the dipole cross sec- 
tion for this pair to scatter off the target at an impact 
parameter b^. The former is well known from QED, 
while the latter represents the dynamics of QCD scat- 
tering at small x. In the "classical" limit without high 
energy evolution effects, the dipole cross section can be 
written as 



d 2 b T 
where T p (by) 



2[1 



e X p(-r 2 F(x,r)T p (b T ))] , (2) 

<-p\"T> ■- ex P( — 2S^j) ^ S ^ ne i m P ac t param- 
eter profile function in the proton, normalized as 
J d 2 h T T p (b T ) = 1 and F is proportional to the DGLAP 
evolved gluon distribution 11] 
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In general, the dipole cross section can range from in 
the r — > color transparency limit to 2, the maximal uni- 
tarity bound. The saturation scale Q s characterizes the 
qualitative change between these regimes; following [l9| 
we shall define here Q s as the solution of 



dcr, 



dip 



d 2 hi 



(x,r 2 = l/Q 2 (x,b T )) = 2(1 -e-V*) (4) 



Note that our definition is completely model independent 
and can be applied to any sensible parametrization of the 
dipole cross section. The definition ^ applied to a Gaus- 
sian dipole cross section gives the saturation scale of the 
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GBW model, but it differs slightly from the convention 
in Ref. [20| , where the saturation criterion was taken as 
-§20 = 2(f-e- 1 / 2 ) at r 2 = 2/Q s 2 . The dipole cross sec- 
tion in Eq. ^ was implemented in the impact parameter 
saturation model (IPsat) [13] where the parameters (j,q, 
C, and Bq (as well as two other parameters characteriz- 
ing the initial condition for the DGLAP evolution) are fit 
to reproduce the HERA data on the inclusive structure 
function F2. 

The form of the IPsat dipole cross section in Eq. ([2]) is 
applicable when leading logarithms in Q 2 dominate over 
leading logarithms in x. At very small x, quantum evolu- 
tion in the CGC describes both the dilute bremsstrahlung 
limit of linear small x evolution as well as nonlinear RG 
evolution at high parton densities [H[ . The essential dy- 
namics of this small x evolution are combined with a more 
realistic 6-dependence in the bCGC model [U, [22]. The 
model is formulated in terms of an explicit x-dependent 
saturation scale that we shall denote by Q' s to distinguish 
it from our model independent definition ([J|. The dipole 
cross section has the form 



2 (7a + ^l n (^)) 



for rQ' s < 2 



= 2 - 2 cxp (-A In 2 (BrQ' s )) for rQ' s > 2 . 

(5) 

The coefficients A and B in the second line of this equa- 
tion can be determined uniquely from the condition that 

d 2b'^ aid its first derivative with respect to rQ' s are con- 
tinuous across rQ' s — 2. Here Y = ln(l/x) is the rapidity, 
while 7 S = 0.628 and K — 9.9 (which quantifies the ge- 
ometric scaling violations in |5])) are obtained from the 
leading logarithmic BFKL dynamics (23j. The impact 
parameter dependence of the proton saturation scale is 
introduced into the bCGC model in the form 

Q' s (x, b) = (^) # [exp (-& 2 /25 CGC )] ^ GeV . (6) 

After choosing A/o = 0.7, the parameters A, xq and Bqgc 
are fit to the data. We will discuss the impact parameter 
dependence of the dipole cross-section and the satura- 
tion scale further in the next section. We must empha- 
size that the saturation scale Q s is conceptually the same 
as Q' s and their numerical values are of the same order, 
but we differentiate between them in order to maintain 
our model independent definition ((4]) and the original 
parametrization of the bCGC model. 

Both the IPsat model and the bCGC model provide 
excellent fits to a wide range of HERA data for x < 
0.01 [H, [H. In Ref. [H[, we discussed the possibility 
that DIS off nuclei can distinguish respectively between 
these "classical CGC" and "quantum CGC" motivated 
models. Our discussion in that paper addressed the A 
dependence of the nuclear saturation scale and fits to the 
available nuclear DIS inclusive data. We also addressed 
elastic scattering of qq dipoles of nuclei but our analy- 
sis was incomplete because a more complete picture of 
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FIG. 1: Kinematics of diffractive DIS. 



diffraction requires that one consider the diffractive scat- 
tering off nuclei of higher Fock states as well — at least 
of the qqg Fock state. This shortcoming is addressed in 
the present paper. We note that there have been several 
discussions of diffractive scattering on nuclei in the liter- 
ature [n in m m, m m m m, hi m, m. we wm 

later compare and contrast our results with those studies 
that overlap with ours. 

The paper is organized as follows. In the next section, 
we will introduce the kinematics for diffractive deep in- 
elastic scattering and the relevant formulae to compute 
diffractive structure functions in the dipole model. We 
will discuss how the qq and qqg dipoles, the dominant 
contributions at present collider energies, contribute to 
the diffractive cross-section. Particular attention will be 
paid to the impact parameter dependence of the cross- 
sections. In Sec. Ill, we will discuss a comparison of 
the different CGC based models to the HERA diffrac- 
tive data. This analysis is extended to nuclei in Sec. IV. 
The "breakup" and "non breakup" events discussed in 
Ref. [l!| for qq dipoles are applied to the qqg dipoles as 
well. We study the nuclear enhancement and suppres- 
sion of diffractive cross-sections in the different models 
and assess their predictive power. In the final section we 
compare our results with those existing in the literature 
and outline future research on this topic. 



II. COMPUTING THE DIFFRACTIVE 
STRUCTURE FUNCTION 

In the dipole picture of diffractive DIS, the virtual pho- 
ton fluctuates into a colorless parton Fock state, which 
interacts elastically with the proton or nucleus. In the 
final state, the latter remains intact, while (for inclu- 
sive diffraction), the parton Fock state fragments into 
hadrons with an invariant mass Mx- In addition to the 
usual DIS invariants x and Q 2 , inclusive diffractive scat- 
tering can be fully characterized by two other invariants 
(5 and t. Here t = -(P- P') 2 , where P» and P' M denote 



3 



the four-momenta of the incoming and outgoing hadron, 
and (3 = Q 2 /{Q 2 + M%), with < (3 < 1. Note that if 
Y = \n(l/x) is the relative rapidity between the electron 
and the hadron, ln(l//3) is the rapidity interval between 
the electron and the hadronic fragments X of the Fock 
state. The rapidity gap between X and the proton or 
nucleus is expressed as Yp = ln(l/xp); therefore, by def- 
inition, Y — ln(l//3) + Yp, or xp = x/ (3. The kinematics 
is illustrated in Fig. [T] 

The parton Fock state could be a qq dipole or higher 
Fock states which involve emission of one or more gluons. 
Naively, these higher Fock states are suppressed by higher 
powers of a s ; however, in the limits of large Q 2 or small 
(3, gluon emissions are accompanied by large logarithms 
in Q 2 and (3 respectively that compensate the powers of 
the coupling. These factors then have to be resummcd 
to obtain the diffractive cross-sections in the appropri- 
ate kinematic domains. The powers of a s In 1//3 can in 
principle be resummed in the dipole picture by using the 
Kovchegov-Levin equation I36H in the mean-field approxi- 
mation, and the approach of ,37j beyond. In practice one 
must require xp < 0.01 for the rapidity gap to be large 
enough to be clearly detectable and to be in the domain 
of validity of the dipole model. Thus in the kinematical 
regime of HERA (and realistic future DIS experiments) 
In 1//3 is not very large and we shall only take into ac- 
count the contribution of the qq and qqg components in 
the present work. 

In diffractive scattering, contrary to inclusive scatter- 
ing, one is really computing the square of an amplitude 
— with interference terms. As a consequence, the square 
of the photon wave function cannot be factorized from 
the cross section. Introducing the auxiliary notation 



d 2 t>7 



drrK n (er)J n (kr) ' P (b T , r, x P ) 
cr b*r 



. (7) 

the qq components of the diffractive structure function, 
can be expressed as 

XfF^ qq (x r , (3, Q 2 ) = y^jg X/ e / J dzz ^ ~ z ) 

x [e 2 (z 2 + (1 - 2 ) 2 )$x + m 2 $ ] (8) 



L,qq( x ViPi 



NcQ 6 V 2 

4tt 3 /3 Y 6f 



1/2 



dzz 3 (l-z) 3 $ 



(9) 



with e 2 = z(l 

f 



z)Q 2 



z(l — z)M\ — m 2 f and 



zo = ^1 — Wl — &m 2 /Mxj /2. In principle the measur- 
able \t\ is bounded by a very small kinematical lower 
bound \t\ > (xpm^) 2 and both upper and lower exper- 
imental limits. Here we have, however, performed the t 
integration from — oo to 0, because this enables us write 
the above expressions in impact parameter space. Since 



the decrease of the structure functions with t is (close to 
an) exponential, this approximation is reasonable. Note 
that we allow for an arbitrary impact parameter depen- 
dence, and thus there is no "diffractive slope" as a sepa- 
rate parameter. 

For the qqg component of the structure function, there 
are essentially two approaches that have been used. The 
one derived in the large Q 2 -limit in Ref. 38] gives 



jnD (GBW) / a 

/ <l 2 b, 



T,qqg 

a s f3 ■ 
8tt 4 
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^ z 

lp (b T ,r T ,a; P 



K 2 (y/zkr)J 2 Wl - zkr) 



(10) 



In this limit, the qqg-system is an adjoint representation 
gg dipole, therefore, as in Ref. [39^ . we use here the dipole 
cross section in the adjoint representation, denoted by 



dad 



ip 



d 2 b 7 



= 2 



111 1 dg dip 

' 2 d 2 b T 



(11) 



This large N c expression, where the t^g-system is treated 
as two fundamental representation qq dipoles, is con- 
sistent with the Balitsky-Kovchegov [4fJ, |4l|, |43] mean 
field treatment of small x evolution. The same is true 
of the j3 — > 0-limit (Eq. (O) below, where the qqg- 
system is also treated as two fundamental representation 
qq dipoles, but this time with different sizes. We will re- 
fer to Eq. (JTUJ) as the "GBW" ggg-componcnt (referring 
to [8j). Note that the form (jTTj) differs from the one used 
in Ref. Q, where the adjoint dipole does not have the 
right saturation limit. 

The other well studied case is the (3 — > limit [H, |43|, 
0, El, HI S3- In this limit (the notation "MS" below 
refers to the authors of 43]) the structure function again 
factorizes into the product of a photon wave function and 
the cross section for the qqg-systera to interact elastically 
(diffractively) with the target: 



i^ ( rw=o,Q 2 ) 



Xpf T 



C F a s Q 2 

47T 4 Q!cm 



d 2 r T / dz 





d 2 b T A(r,x P ,b T ), (12) 



with 



A(r, xp,b T ) = J d 2 r T '- 



r T 



yt' 2 {yt — tt' 



jV(r T ') 



I 2 



+ jV(tt - r T ') - AA(r T ) - Af(r T r )Af(r T - v T ') 



(13) 
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FIG. 2: /3-dependence of the different contributions to the 
proton diffractive structure function at Q 2 — 5 GeV 2 and 
x w = 10~ 3 . 



and Af(r, xp,h T ) — (r, xp, h T )/2. The coupling con- 
stant a s in Eqs. (fT0|) and (fT2|) is treated here as a con- 
stant free parameter (independent of the DGLAP evolu- 
tion momentum scale in the IPsat model). A more thor- 
ough study of the running coupling effects in this problem 
is an interesting question that is out of the scope of the 
present work. 

Depending on the mass of the diffractive system Mx 
or, equivalently, /3, the diffractive structure function is 
dominated by either the qq or qqg Fock states. (Sec 
Ref. [HI for a general argument of the j3 dependence.) 
Specifically, in the limit (3 — > 1, the dominant component, 
^Lqq 1S a longitudinally polarized qq system. At interme- 
diate /3 ~ 0.5 the dominant component is a transversally 
polarized qq denoted here by ^. In the limit (3 — > 
the invariant mass of the diffractive system is large, and 
this large phase space is filled by radiation of additional 
gluons, each of them being suppressed by a s . This struc- 
ture is illustrated in Fig. O 

In Ref. 39], it was shown that the (3 — limit 
of Eq. (fT0|) , at Q 2 — > oo, approaches the result from 
Eq. (fT2| . This therefore suggests the following interpo- 
lation formula between the two limits (39|: 



(GBW) 

X ^ F T,qqq 



(x P ,/3 = 0,Q 2 ) 



(14) 



We also note the work in Refs. (3CJ, [32|, [4!| , which uses 
yet another prescription whose relation to ours is not very 
transparent. 



III. IMPACT PARAMETER DEPENDENCE 

We shall now discuss the ^-dependence of the dipole 
cross-sections (see also e.g. Refs. [5(1 [U [Hj]). Several 




FIG. 3: 6-dependence of the the inclusive structure func- 
tion and different contributions to the diffractive structure 
function at Q 2 = 1 GeV 2 (a) and Q 2 = 100 GeV 2 (b) for 
x — 10~ 3 (inclusive) and xp = 10 -3 (diffractive). In plot 
(b) (Q 2 = 100 GeV 2 ) the 6-dependence of the inclusive cross 
section and gggr-components are indistinguishable. 



works on the subject (for example, Refs. [1, [U HH) 
assume, explicitly or implicitly, a factorizable hx depen- 
dence 

<Wdip -(b T) r T , : r) = 27v'(bT,r T ,x)=2T p (bT)iV(r T , : r). 



d 2 b T 

When considering diffractive scattering on protons, this 
is consistent with the exponential t dependence observed 
in experiments, and in fact implies that Tp(b^) is Gaus- 
sian. 

In the IPsat model, in contrast to the factorization 
of the 6-dependence in Eq. (fl~5|) . the dependence of the 
dipole cross-section on impact parameter is as in Eq. © , 
which equivalent to an impact parameter dependence of 
the saturation scale Q s 2 oc T p (br). In the IPsat model 
the impact parameter profile of the proton saturation 
scale is chosen to have the form 



T P (h T ) 



1 



2ttB g 



(16) 



which is normalized to unity. In the large <5 2 -limit the 
cross section is dominated by small dipoles and one can 
expand the exponential of Eq. |(2]), and the dipole cross 
section becomes proportional to TLQby). This corre- 
sponds to (b 2 ) = 2Bq, which can be interpreted as the 
average square gluonic radius in the proton. 
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We shall now justify in the framework of the BK equa- 
tion why we prefer (nj [2(| [HI to incorporate the impact 
parameter dependence in the saturation scale Q s instead 
of as a multiplicative factor. (As we shall see, this dif- 
ference is more significant for diffractive as opposed to 
inclusive observables.) The fully impact parameter de- 
pendent BK equation for the dipole amplitude is 



a s N c 
2tt 2 



d zt- 



xt - zt) 2 (vt - zt) 2 



(17) 



where t>T = (xj- + yr)/2 and rr = xt — yr- This equa- 
tion was studied numerically in Ref. [53| with the conclu- 
sions, a) that the t>T and yt dependence of the amplitude 
does not factorize, and b) one obtains unphysically large 
power law tails in by. Physically this is due to the fact 
that the BK equation does not include confinement ef- 
fects that would cut off these tails. An approximation 
that removes the growth of the power law tail with in- 
creasing y — log(l/x) is based on the argument that the 
dipole sizes r' ~ 1/Q S (which dominate the integration in 
(jTTJ) ) are parametrically smaller than the typical scale R p 
for the variation of the amplitude with impact parameter 
b. One can thus approximate (note that the same approx- 
imation was already implied when writing Eq. (| 13() ) 



diffractive final states 1 . At smaller Q 2 (i.e. larger r) 
the impact parameter profile in the IPsat model is not 
exactly Gaussian, and diffractive peaks appear [20| at 
large —t. This is, however, in a region which is most 
likely unobservable because of the exponential suppres- 
sion of the diffractive cross section. We emphasize (see 
also Ref. [l!|) that specifying the full bT-dependence of 
the dipole cross section simultaneously determines both 
the normalization ("oo") °f the cross section and the 
diffractive slope B — dlna D / di|t = o- 

In Fig. [21 we plot the inclusive and diffractive struc- 
ture functions as a function of 6, normalized to the same 
quantity at b — 0. As has been noted previously, the in- 
clusive distribution does not change significantly as Q 2 is 
varied. The qqg dipole contribution, at low Q 2 , is signif- 
icantly broader in b than the inclusive distribution. This 
should be interpreted as its relative suppression at small 
Q 2 and small b, where the dipole cross section is close to 
the black disk limit and the qqg contribution vanishes, as 
is easily seen from Eq. (flU)) . 

Inclusive diffraction off nuclei is not very sensitive to 
the choice of impact parameter profiles in the proton. 
However, there are significant issues related to the impact 
parameter distribution of nuclei that are very relevant 
for diffraction, as we will see later. But before moving to 
scattering off nuclei we will first discuss the comparison 
of the IPsat and bCGC models to the HERA data on 
inclusive diffraction in ep scattering. 



dV 



2tt 2 J 1 r T ' 2 (r T -r T ') 2 
N (b T , r T ') + J\f (b T , r T - r T ') ~ A/" (b T , r T ) - 

AT(b r ,r T ')A/'(b T ,r T - r T ') 



(18) 



It is obvious that if the 6-independent equation has a 
scaling solution that can be expressed in terms of a sat- 
uration scale Qs, then replacing it with an impact pa- 
rameter dependent Q s (i>T) gives a solution of Eq. (fTSl) . 
This is effectively the approach used in Refs. [3ll[5^.l53|. 
The factorized ansatz Eq. (TT5]) is not (unless the profile 
is a 8 function) a solution of even this approximate equa- 
tion. A factorized Gaussian profile for the proton dipole 
cross section, for example, does not approach the correct 
unitarity limit for 6^0. This is the main argument for 
including the impact parameter dependence in the satu- 
ration scale, not as a factorizable prefactor of the dipole 
cross section. 

Because b is the Fourier conjugate variable of A, where 
t = — A 2 is the momentum transfer squared between the 
incoming and outgoing proton, a Gaussian dependence of 
the impact parameter profile corresponds to an exponen- 
tially decreasing dependence of the dipole cross-section 
on t. This exponential behavior is widely observed for 



IV. COMPARISON TO HERA DATA 

An extensive comparison of the b-CGC and IPsat 
dipole models to the HERA data on exclusive vector me- 
son production was performed in Ref. [12] . We compare 
our calculation to the HERA results on diffractive struc- 
ture functions for xp < 0.01, measured both using the 
rapidity gap method (ZEUS FPC [tf| and HI LRG |3) 
and by measuring the recoil proton (ZEUS LPS [58| and 
HI FPS [13). Because the FPC and LRG data include 
events in which the proton has broken up, the cross- 
sections measured for the process ep — > eXY are larger 
than the one measured for the process ep eXp. We 
scale down this data by a constant factor to correct for 
the proton dissociation contribution; the ZEUS FPC data 
by a factor of 1.45 and the HI LRG data by 1.23. These 
factors are different due to the different cuts on My , the 
mass of the proton system. Note that it is the FPS- 
LPS data that correspond to our definition of diffractive 
events and to our formulae, as the proton should escape 
the collision intact. These correction factors were deter- 
mined in the region < 0.7 where data exists from both 
methods and then assumed to be independent of (3; the 



In contrast, the profile T(by) = 
distribution [H] ~ 4| Jl(^/=t.R)| 2 /(- 



26(R - b), 
-tR 2 ) «14 



leads to the t- 
tR 2 /A + 0{t 2 ). 
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FIG. 4: Comparison of the IPsat and bCGC fits to HERA data on diffractive structure functions. 



largest contributions to the \ 2 of our analysis come from 
the large (3 region where this assumption has not been 
tested experimentally. All the x 2 -vahies quoted below 
are calculated per degree of freedom using statistical and 
systematical errors added in quadrature. The total num- 



ber of experimental points in the four data sets is 343, of 
which 76 are from the FPS-LPS data. 

The experimental data are presented in terms of the 
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reduced cross section a®' 3 (l3,xp,Q 2 ) 

i-y + y ) K [p,%¥,Q ) 



d/3 dQ 2 



PQ 4 



a 



D,3 j^D,3 



2 - 2y 23,3 

-Cr 



2-2y + y 2 



2 L 



(19) 
(20) 



with y = Q 2 /(sx) where y / s = 318 GeV is the total energy 
in the e—p collision. 

With the IPsat cross sections, for the combined dataset 
from ZEUS and HI data both with and without identified 
protons we get \ 2 — 1-3 for a s = 0.14 in Eqs. (|10p and 
(Ti"2"]) . For the bCGC cross section the fit is equally good 
and works with a larger value of the coupling: \ 2 = 1.3 
with a s = 0.22. These are the values of a s that we shall 
use for the respective models to evaluate nuclear diffrac- 
tive structure functions in the next section. For the IPsat 
model the largest contribution to the \ 2 comes from the 
rapidity gap method data at large f3. The fit to only 
the LPS (x 2 = 0.5 IPsat) and FPS ( X 2 = 0.8) is much 
better. Considering just the LPS also accommodates a 
larger value of a s — 0.21 with still \ 2 < 1. The fit for 
bCGC is more even among the datasets, but also there 
the HI rapidity gap data has a larger \ 2 — 1-9 than the 
other data sets. Our fit of these two models to a combi- 
nation of the HERA datasets is presented in Fig. [U 

The fit to HERA data is better with a smaller a s than 
in Ref. (3{|. Given the ^-dependence described previ- 
ously this is to be expected. The factorized ^-dependence 
(Eq. ([To]) ) used in earlier calculations of the diffractive 
structure function such as Refs. [1, H^] forces the qqg- 
component to have the same impact parameter depen- 
dence as the (yg-component. In a more realistic descrip- 
tion, as suggested by Fig. [31 the qqg component is sen- 
sitive to larger impact parameters and is thus larger; in 
order to fit the same data this must be compensated by 
multiplying it with a smaller factor of a s . 



V. THE NUCLEAR DIFFRACTIVE 
STRUCTURE FUNCTION 

A straightforward generalization of the dipole formal- 
ism to nuclei is to introduce the coordinates of the indi- 
vidual nucleons {b^}. In the IPsat model the replace- 



ment Tp(by) 
section 



J2i=i T p (br — bri) gives the dipole cross 



d^dip 
d 2 b T 



= 2 



1 _ e -r 2 F(x,r) Ef =1 T p (br-b-z 



(21) 



where F is defined in Eq. (J3J). The positions of the nucle- 
ons {b^} are distributed according to the Woods-Saxon 
distribution Ty^b^i). We denote the average of an ob- 
servable O over {b Ti } by 



A 

j n d 2 b Tl T A (b Tl )o({b Tl }) 



(22) 
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FIG. 5: The ratio F$ A X /(AFf p x ) in the IPsat model for dif- 
ferent components of the diffractive structure function (a; = 
T, L : qqg) as a function of (3 for gold at Q 2 = 5 GeV 2 and 
xp — 10 -3 without nuclear breakup. 
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FIG. 6: Same as Fig. [5] but including nuclear breakup. 



Here we have introduced the Woods-Saxon thickness 
function 



T A {b T ) = [ dz- 
J 1 



C 



exp 



(v/b 



T z + z 2 -R A )/d 



(23) 



which is normalized to unity J d 2 hTT A (hT) = L The 
nuclear radius R A and surface diffuseness d are measured 
from the electric charge distribution, their values can be 
found in Ref. [60| . The average differential dipole cross 
section is well approximated by [20| 



d 2 b T 



1 T A (h T ) p 
1 2 di P 



(24) 



where, for large A, the expression in parenthesis can be 
replaced by form |30j 



exp 



AT A {b T ) 



P 

^dip 



(25) 
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FIG. 7: The ratio F^ A /{AFQ as a function of for Ca, Sn 
and Au nuclei for Q 2 — 5 GeV 2 and — 10 -3 . Results are 
for the "non breakup" case in the IPsat model (thick lines) 
and the bCGC model (thin lines). 



All parameters of the model come from either fits of the 
model to ep-data or from the Woods-Saxon distribution; 
no additional parameters are introduced for eA collisions. 

The Glauber form (|24f has a straightforward in- 
terpretation as the dipole scattering independently off 
the different nucleons. To see this explicitly denote 
-^^(r T ,b T ) = 2(1 - S(rT,br)), where the S-matrix 
element S(tt, b-r) is the amplitude for the dipole to 
not interact (elastically; the relation to the inclusive 
cross section is via the optical theorem) with the tar- 
get. The S- matrix element for scattering off a nucleus 
is then given by S A (r T ,h T ) = l\f =1 S p (r T ,b T - b Ti ) 
which, for the IPsat model, turns out to be equivalent 
to T p (b T ) -> J2t=i T pi h T - b Ti ). Note that in the form 
(|24p there is no leading twist shadowing, i.e. in the large 



Q or small r limit cr 



dip 



because in this limit 



crj ip ~ r is small and one can expand the exponential. 

The situation for the bCGC model is much more com- 
plicated, since the replacement Tpfox) — > J2i=i T p (br — 
hxi) into the definition of the bCGC saturation scale 
© does not lead to the Glauber form (|24p. One could 
see this as a consequence of the "noncommutativity" of 
nuclear effects and high energy evolution: even if one as- 
sumes that for a particular x and vt a dipole interacts 
independently with the nucleons in a nucleus, this will 
not necessarily be the case for other rapidities and dipole 
sizes because the evolution sums up nonlinear interac- 
tions between the nucleons. Since it is not completely 
obvious how to introduce a nuclear dependence directly 
into the bCGC parametrization for the dipole cross sec- 
tion we will in this work use (|24p for the bCGC model as 
well. A comparison of high energy evolution for protons 
and nuclei would be out of the scope of this work, see 
however Refs. [6l|, [62[ . 

In Ref. [l!|, we showed that the nuclear dipole cross- 
sections obtained in this manner gave a good (parameter 
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FIG. 8: The ratios F% A X /{AI% X ) at x r = 10 -3 for different 
components of the diffractive structure function plotted as a 
function of Q 2 . The components are evaluated where they 
are dominant: at /3 — 0.1 for qqg, /3 = 0.5 for T and /3 = 0.9 
for L. Results are in the IPsat model for both "breakup" and 
"no breakup" cases, (a) Ca nuclei, (b) Au nuclei. 



free) agreement with the x and Q 2 dependence of the 
NMC inclusive structure function data [63j, [64( at small 
x. However, at the level of the accuracy of the data, it 
was not possible to distinguish between the IPsat and 
b-CGC models for the inclusive cross-section. We will 
now consider nuclear diffractive (qq and qqg) structure 
functions in the two dipole models. This is obtained by 
substituting the nuclear dipole cross-section (Eq. (|24|) ) in 
Eqs. ©, GUI and 02). 

It is very easy to break up a nucleus with a rela- 
tively small momentum transfer \t\ > |f^; n |. However, 
for |tJn in | ^ |i| ^ Kminli where i^ ln is the minimum mo- 
mentum transfer required to break up the proton, one 
can still have a nuclear diffractive event with a rapid- 
ity gap. For \t\ > \t^ in \, the "lumpiness" of the nu- 
cleus shows up as a proton-like tail ~ exp{CtRp)} of 
the ^-distribution. In our formalism, if one requires that 
the nucleus stays completely intact, the average (-) N in 
Eq. (12"2")) must be performed at the amplitude level; the 
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results are therefore proportional to 



One 



finds that in this case da D / dt falls off very rapidly as 
~ exp{CtR A }. We refer to these (known as coherent 
diffraction) as "non breakup" events. Measuring the in- 
tact recoil nucleus at such a small t experimentally at 
a future electron ion collider [65| is challenging. Also 
including events where the nucleus breaks up into color 
neutral constituents without filling the rapidity gap be- 
tween the qq dipole and the nuclear fragmentation re- 
gion [incoherent diffraction) corresponds to performing 
over the cross section (66|; in this case, 

) instead. The 



the average 

one performs the average 



d 2 b T 



difference between the two averaging procedures can be 
significant with increasing values of t as shown explicitly 
in Ref. (l9j |. The "breakup" figures in our plots are a sum 
of these incoherent and coherent events. 

In Figs. and [5] we show the ratios of different com- 
ponents of the gold diffractive structure function to the 
proton one as a function of (3. The qq-components of 
the F^ A are enhanced compared to A times the proton 
diffractive structure functions. This is to be expected, 
because of the fact that in a gold nucleus the dipole cross 
section is, on average over the transverse area, closer to 
the unitarity limit than the proton - it is "blacker" . The 
elastic scattering probability of a qq dipole is maximal in 
the "black" limit and the approach to it is quicker in a 
large nucleus. The qqg component, on the other hand, 
is suppressed for nuclei compared to the proton. This is 
due to the fact that in a nucleus the scattering ampli- 
tude is closer to the unitarity limit when the qqg compo- 
nent vanishes, as can be seen e.g. from Eq. (fT5|) . This 
leads to a nuclear suppression of the diffractive structure 
function in the small (3 region, where the qqg component 
dominates. The net result of the different contributions 
is that F^ A , for a large range in (3, is close to AF^. 

In Fig. [7J we plot the total ratio as a function of (3 for 
different nuclei in the "non breakup" case. As expected 
from our prior discussion, one sees a strong enhancement 
with A for larger (3 and likewise, a stronger suppression 
with A at very small values of (3. A comparison of the 
"breakup" versus "non breakup" cross-sections can be 
seen in Fig. [5] for the ratio of diffractive cross-sections as 
a function of Q 2 . The results in Fig. [5] for the ratio of 
diffractive structure functions indicate that the diffrac- 
tive cross-section in nuclei decrease more slowly for large 
Q 2 than in the proton. 

In Figs.l9"l[T0"land[TTlwe show the xp dependence of the 
nuclear modifications. For a fixed Q 2 — 5 GeV 2 the nu- 
clear enhancement of the qq components becomes smaller 
at smaller xp. This can easily be understood as an anal- 
ogous effect to the Q 2 -dependence plotted in Fig. [5] in- 
creasing xp for a fixed Q 2 increases Q 2 /Q s 2 (xp) and has 
the same effect as increasing Q 2 at a fixed xp. In Fig. [TU] 
we compare the xp- dependence in the IPsat and bCGC 
models. As was already observed in [l9l |. the experimen- 
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tal signature for the different evolution dynamics in the 
models in eA-scattering is mainly in the different xp de- 
pendence in the nuclear modification factor. The result 
presented in Fig. [10] confirm this, although the effect is 
perhaps smaller than in the shadowing of the inclusive 
cross section. In Fig.[TT]we compare, in the IPsat model, 
the nuclear modifications to F® q q at f3 = 0.5 for Ca, 
Sn and Au nuclei. One can see that the "non breakup" 
curves are much more sensitive to the nuclear species 
than the "breakup" ones. 

In Fig. [T21 the dependence of the longitudinal and 
transverse components of the diffractive structure func- 
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straightforward to give a very precise general estimate 
for <7Tj/<7 to t that would be observed in a generic high 
energy eA collider. A general order of magnitude ar- 
gument would be as follows. For moderate values of Q 2 
and large nuclei we expect a nuclear shadowin g of the 
inclusive structure function by a factor ~ 0.8 [lj|. A 
typical nuclear enhancement of diffraction (at moderate 
values of [3 > 0.2) is a factor of ~ 1.2 (see e.g. Fig. [7]). 
Combining these we expect au/o'tot to be increased by 
a factor of 1.2/0.8 = 1.5 compared to the proton. Thus 
from a typical ep fraction of 15% we expect ovjotot to 
go up to 20% - 25% at an eA collider. 



VI. COMPARISONS AND CONCLUSIONS 



FIG. 11: Fg q9iA /AFg q g tP at Q 2 = 5 GeV 2 and = 0.5, as 
a function of xw for Au, Sn and Ca nuclei in the IPsat model. 
Both breakup and non-breakup cases are shown. 
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tion on nuclear size is shown for the "breakup" and "non 
breakup" cases. In the "breakup" case, one sees a very 
weak A dependence In the coherent "non breakup" case, 
one first notes that the diffractive structure function first 
decreases up to atomic numbers A ~ 10, before begin- 
ning to rise. As noted in Ref. [19], this is due to the 
typical scattering amplitude for small nuclei actually be- 
ing smaller than for a proton because of the diluteness 
of the nucleus. This leads to a suppression of coherent 
diffraction. The "breakup" case, on the other hand, can 
only be enhanced in nuclei. For gold nuclei, the cross 
sections in the "non breakup" case are about 15% lower 
than in the "breakup" case. 

Because of the different nuclear modifications in in- 
clusive and diffractive scattering, the fraction of diffrac- 
tive events in an experiment depends on the detailed 
kinematics and experimental coverage. Thus it is not 



Finally we shall briefly compare our results to some 
other recent literature on diffraction in eA scattering. 
Several aspects of the treatment of nuclear diffraction 
here were discussed previously in the works of Nikolaev, 
Zakharov and Zoller [2(| and of Frankfurt and Strik- 
man [27| — these works may be consulted for earlier 
works in the literature as well. With regard to these 
early works, a key difference is the use of the explicit 
form of the qqg contributions to the diffractive struc- 
ture functions given in Eq. (TTUJ) and (TT2"|) . Further, we 
have focused on the predictions of Color Glass Conden- 
sate based models, that were previously constrained from 
fits to HERA data. In this paper, we have made explicit 
fits of these models to the HERA diffractive structure 
function data as well for the first time. 

Turning to relatively more recent works in the litera- 
ture, Gotsman, Levin, Lublinsky, Maor and Tuchin [30| 
use a parametrization of the impact parameter depen- 
dence that is similar to ours to compute the fraction of 
diffractive events in the total cross section. They see a 
modest enhancement compared to ep scattering in agree- 
ment with our results. As previously mentioned, they use 
a different method to calculate the qqg contribution but 
their results (see for example, Fig. 7 in Ref. [1^) seem to 
point towards its relative suppression in nuclei, in qual- 
itative agreement with ours. On the other hand the en- 
ergy dependence seems much stronger than our results 
indicate. Since the results of Ref. [3fJ are not presented 
in the form of diffractive structure functions (in particu- 
lar the /3-dependence is not calculated) a more detailed 
comparison is difficult. In Ref. [32j Levin and Lublinsky 
start from a Glauber-like parametrization very similar 
to our Eq. I|25p. but end up with a saturation scale de- 
pending on the nuclear size as Q s 2 ~ A 6 for moderately 
small values of x ~ 10 -3 , which naturally leads to a much 
stronger A-dependence of (To/atot than our results. 

In the work of Frankfurt, Guzey and Strikman [28l . 
l29l the nuclear diffractive structure functions are 
modelled starting from leading twist diffractive parton 
distributions in the proton. Although the terminology 
and theoretical framework are somewhat different, some 
comparisons can still be made. They find that F^ A /A is 
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suppressed compared to F^, as a function of Q 2 , both 
at small (3 = 0.1 (in agreement with our result) and in 
the transverse qq dominated region (3 — 0.5 (which dis- 
agrees with our findings); this can be seen explicitly by 
comparing Fig. 5 in [67| and our Fig. [8] 

Kugeratski, Goncalves and Navarra [H, [H, [6§| com- 
pute nuclear diffractive structure functions using the 
GBW framework Q , albeit with the IIM [2l[ dipole cross 
section. They extend the calculation to nuclei using a 
simple A 1 / 3 -scaling of the saturation scale, the trans- 
verse area a and the diffractive slope Bjj; the latter 
two are treated as independent parameters. They ob- 
serve (see Figs. 4 and 6 in (35|) a similar pattern in the 
/3-dependence as we do, namely that the nuclear F% is 
relatively suppressed at small (3. However, from their 
results one also deduces (note that the diffractive struc- 
ture functions in [3|| are normalized by A 4 / 3 which is 
not explicit in the text [z3]) that F^ A /AF^ p is of the or- 
der of 1/4, which is much smaller than our result ~ 1 
(see Fig. [7]). We conjecture that these differences arise 
because, unlike in our discussion, an and Bp are treated 
as independent parameters in Ref. [35j, with values that 
are not necessarily consistent with each other. 

Measuring diffractive events in a high energy eA col- 
lider would yield significant insight into the physics of 
high parton densities. The experimental measurements 
of inclusive diffraction are encoded in diffractive struc- 



ture functions that depend on the kinematical variables 
of the process. We have argued that it is possible to make 
a large number of detailed theory predictions for nuclear 
diffractive structure functions. Comparing these to ex- 
perimental measurements would provide a powerful han- 
dle to understanding the behavior of QCD at high energy. 
In this paper we have presented some of these theory pre- 
dictions, within two particular models for the high energy 
wavefunction. One (IPsat) is constructed to combine the 
correct DGLAP limit at large Q 2 with unitarization at 
small Q 2 . The other (bCGC) is a parametrization of the 
leading order BK dynamics. We have also emphasized 
the importance of including in the calculations a realistic 
description of the impact parameter dependence in the 
nucleus, without which an actual comparison to experi- 
mental data is not very useful. 
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